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Abstract We have studied two types of meshwork models by using the 
canonical Monte Carlo simulation technique. The first meshwork model 
has elastic junctions, which are composed of vertices, bonds, and trian- 
gles, while the second model has rigid junctions, which are hexagonal (or 
pentagonal) rigid plates. Two-dimensional elasticity is assumed only at the 
elastic junctions in the first model, and no two-dimensional bending elastic- 
ity is assumed in the second model. Both of the meshworks are of spherical 
topology. We find that both models undergo a first-order collapsing tran- 
sition between the smooth spherical phase and the collapsed phase. The 
Hausdorff dimension of the smooth phase is ~ 2 in both models as ex- 
pected. It is also found that ~ 2 in the collapsed phase of the second 
model, and that H is relatively larger than 2 in the collapsed phase of the 
first model, but it remains in the physical bound, i.e., H <3. Moreover, the 
first model undergoes a discontinuous surface fluctuation transition at the 
same transition point as that of the collapsing transition, while the second 
model undergoes a continuous transition of surface fluctuation. This indi- 
cates that the phase structure of the meshwork model is weakly dependent 
on the elasticity at the junctions. 
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1 Introduction 

Two-dimensional curvature model of Helfrich and Polakov as well as that 
of Nambu and Goto has been extensively studied as a model of membranes 
and of strings from the viewpoints of two-dimensional differential geom- 
etry and statistical mechanics [Il[2ll3l|4] . Because of their two-dimensional 
nature, these models have a variety of surface shapes and the correspond- 
ing shape transformations or phase transitions, which can be considered to 
represent the complexity of real physical membranes |5l[6ll7] . 

Among the interesting topics on those models, the surface crumpling 
phenomenon is an old topic that has long been studied from several per- 
spectives [StiQ ulO. . The phase transition of such crumpling phenomena is 
itself an interesting topic in biological and artificial membranes. Experi- 
mental investigations show that such phenomena can be seen in an arti- 
ficial membrane In the string model context the path integration of 
the model describes the sum over surfaces in p]2l[T3] . and therefore it 
seems that the summation technology changes depending on whether the 
surfaces are smooth or crumpled. The crumpling phenomena can also be 
seen in thin sheets. The universal structure was found in the formations 
of singularity in ridges and cones on those crumpled sheets [MlllSj . More- 
over, if the transition is of second-order, the phenomena can be linked to a 
universal model for two-dimensional systems [16^ . 

Thus, extensive numerical studies on the phase transition have been 
made on triangulated surfaces [T7 ( fT8 l fT9 p 20 1 [2T] . including self-avoiding ones 
[22ll23ll24] . Recent numerical studies show that the conventional surface 
model has a first-order surface fluctuation transition on triangulated fixed- 
connectivity spheres pSll^MTTj . and moreover, that the surface collapsing 
phenomena occur at the same transition point; the collapsing transition is 
considered to be a first-order transition. 

On the other hand, the conventional homogeneous surface model men- 
tioned above can be extended to an inhomogeneous one by including a cy- 
toskeletal structure PSII^IHD] . In fact, the phase structure of skeleton mod- 
els is partly understood [29j[30| . The numerical results show that the phase 
structure of the surface fluctuation phenomenon in the fixed-connectivity 
conventional surface model remains almost unchanged even when the com- 
partmentalized structure was introduced |29| . On the contrary, in a dy- 
namically triangulated fluid surface model the phase structure considerably 
changes if the free diffusion of vertices is confined inside the compartments, 
which are introduced as an inhomogeneous structure [30j . 

The compartmentalized models are those defined on a triangulated lat- 
tice with a cytoskeletal structure, which is a sublattice. Thus, the compart- 
mentalized model is defined by using the lattice (= the surface) and the 
sublattice (= the cytoskeleton) , because the one-dimensional bending en- 
ergy is defined on the sublattice and the Gaussian bond potential is defined 
all over the lattice including the sublattice. 
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Therefore, it is natural to ask whether the surface shape is maintained 
only by the cytoskeletal structure. The problem which we consider is whether 
it is possible to eliminate the surface from the compartmentalized surface 
or not. If this is possible, then we are interested in whether the resulting 
model is well defined or not in the sense that the two-dimensional surface 
structure remains unchanged. Moreover, it is also interesting to see whether 
the phase structure remains unchanged when the surface is eliminated. 

Following to these considerations, we studied a meshwork model in [ST] 
and reported some preliminary results on the phase structure. In this pa- 
per, we study two types of meshwork models including the one in [31 on 
relatively large sized surfaces. The first model in this paper is character- 
ized by elastic junctions and is identical to the model in [Sl^, while the 
second model is characterized by rigid junctions, which are hexagonal (or 
pentagonal) rigid plates. 

We will see a first-order transition of surface-collapsing phenomena in 
both models. With respect to surface-fluctuations, the first model under- 
goes a first-order transition, while the second model a second-order transi- 
tion. Thus, our results show that the phase structure of meshwork model 
is weakly dependent on the elasticity at the junctions. The Hausdorff di- 
mension H of the surface can be defined by using the mean square size 
and the total number of vertices N such that ^ N"^/^ . It will also be 
shown that 7?~2 in the smooth phase, and that H remains in the physical 
bound, i.e., H <3, even in the collapsed phase in both models. 

Thus, the first model in this paper has almost the same phase struc- 
ture as that of the fixed-connectivity conventional surface model [26]. On 
the contrary, the phase structure of the second model is slightly different 
from the fixed-connectivity conventional model because of the continuous 
transition of surface fluctuation. 



2 Models 

The meshwork, as mentioned in the introduction, is constructed as com- 
partments on the triangulated spherical surfaces by eliminating vertices, 
bonds, and triangles inside the compartments. Thus, the triangulated sur- 
faces for constructing the meshwork are identical to the lattices for the 
surface models in [29,30J. Therefore, the size of meshwork can be charac- 
terized by the expression similar to the one for those compartmentalized 
surface models in [29ll30] . 

Let N be the total number of vertices, Ns the total number of vertices 
on the chains, and Nj the total number of junctions. Thus, the meshwork 
size can be denoted by {N, Ns , Nj , L) , where L is the total number of bonds 
in a chain between the junctions; L — 1 is the total number of vertices on 
the chain. 
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(a) A meshwork with elastic (b) A meshwork with rigid 

junctions (model 1) junctions (model 2) 

Fig. 1 (Color online) (a)A meshwork with elastic junctions of size 
{N,Ns,Nj,L) = (2562, 1440, 162,4), where iV is the total number of vertices, 
Ns is the total number of vertices on the chains, Nj is the total number of junc- 
tions, and L is the length of chains between junctions, and (b) a meshwork with 
rigid junctions of size (iV, iVg, iVj, L) = (1602, 1440, 162, 4), where N = Ns + Nj. 
The surfaces are drawn in order to visualize the meshwork clearly. 

Figure [^a) is a meshwork of size {N, Ns,Nj, L) = (2562, 1440, 162, 4) 
of the first model (denoted by model 1), where N includes the vertices in 
the junctions, and Nj includes 12 pentagonal junctions. The small dots on 
the chains denote the vertices. The surface of the sphere is shown in the 
snapshot in order to clarify the meshwork. Figure [H^b) shows a meshwork 
of size (TV, Ns, Nj, L) = (1602, 1440, 162, 4) for the second model (denoted 
by model 2), where N is the total number of vertices, Ns, Nj, and L are 
identical to those of model 1. The junctions are counted as vertices in N 
and hence, N is given also by N = Ns + Nj in model 2. 

The meshwork is constructed as follows: Every edge of the icosahedron 
can be divided into £ pieces of uniform length, and then we have a trian- 
gulated surface of size A'o — IQi^ -I- 2 (= the total number of vertices on 
the surface) . The compartmentalized structures are obtained by dividing £ 
further into m pieces (m = 1, 2, • • •). Thus, we have the chains of uniform 
length L = {£/m) — 2 when m divides £. The reason for the subtraction 
—2 is because of the junctions at the two end points of the chain. On the 
meshworks in Figs. HJa) and Wi^), ^ and L are given by £ = 24 and £ = 4, 
and therefore m = 4. 

By using two integers £ and m, we have Nj = lOm^ -I- 2 and Ns = 
30m{£ — 3m) in both models, and therefore A^ = 30to^ — 20m^-|-2 in model 
1 and N = 3Qm£ — 80m^ + 2 in model 2. Since the junctions are considered 
as vertices of the sublattice, we have the expression of Nj = lOrn^ + 2. The 
total number of bonds in the sublattice is 3Nj—6, and each bond contains 
L—1 vertices, then we have A^s = (3iVj— 6)(L — 1). By using L={£/m)—2, we 
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have the above expression of Ns- Nj junctions in model 1 contain 7Nj~12 
vertices and therefore, we have N — Ns + 7 Nj — 12 in model 1, whereas Nj 
junctions in model 2 contain Nj vertices and therefore, we have N = Ns+Nj 
in model 2. 

The total number of the compartments depends on the size N, and in 
fact, it increases with increasing N. However, the chain length L can be 
chosen to be constant and independent of N. We fix the chain length L to 

L = 4, (1) 

which corresponds to n= 10, which is the total number of vertices inside a 
compartment of the surface model in Rcf . [30 . 

Figure [2l[a) shows a hexagonal elastic junction of model 1. The unit 
normal vector is defined on the triangle i at the junctions of model 1, 
and the angle is defined not only on the vertices of the chains but also 
on the corners of the junctions in both models. Figure EJb) shows a rigid 
junction of model 2. 




(a) an elastic junction (model 1) (b) a rigid junction (model 2) 

Fig. 2 (a) A hexagonal elastic junction of model 1, and (b) a hexagonal rigid 
junction of model 2. The unit normal vector in (a) is defined on the triangle 
i in the hexagon, and the angle 9(ij) in S2 is defined not only on the vertices of 
the chains but also on the corners of the junctions in (a) and (b) . 

The Hamiltonian S{X) of model 1 is given by the linear combination of 
the Gaussian bond potential Si, the one-dimensional bending energy 5*2, 
and the two-dimensional bending energy Sj. Si is defined not only on the 
chains but also on the junctions, S2 is defined only on the chains, and Sj 
is defined only at the junctions: 

S{X) = Si + bS2 + bjSj, 5i =^(X,-X,)', 

iv) 

•^2 = [1-^°^%^)] ' Sj^Y.il~n,-n,), (model 1), (2) 
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where Xi{G R^) denotes the three-dimensional position of the vertex i. 
'^(^ij^ in Si is the sum over the bond (ij) connecting the vertices i and j 
on the chains and on the junctions, and i^ '^'2 is the sum over bonds 

i and j, which contain not only bonds in the chains but also bonds that 
connect the center and the corners of the junctions. in Sj is the sum 

over triangles i and j, which share the center of the junction as the common 
vertex. The symbol in S2 is the angle between the bonds i and j, and 
in Sj is the unit normal vector of the triangle i at the junctions as shown 
in Figilja). 

The coefficient b is the one-dimensional bending rigidity {— coefficient 
of one-dimensional bending energy), which will be varied in order to see 
the phase structure, and bj is the two-dimensional bending rigidity at the 
junctions. In this paper, bj is fixed to 

bj = 5 (3) 

so that the junctions are sufhciently smooth. The value bj = 5 is relatively 
larger than the first-order transition point be — 0.77 in the fixed-connectivity 
surface model p51l27j . Therefore, the hexagonal or pentagonal junctions are 
almost flat even when the meshwork is in the crumpled phase at sufficiently 
small b. 

Model 2 is defined on the meshwork with rigid junctions, such as the 
one shown in Fig[TJb).The Hamiltonian is given by the linear combination 
of the Gaussian bond potential 5*1 and the one-dimensional bending energy 
S2 such that 

S{X)^Si+bS2, 

5*2 ^ [1 - COS ^(y)] , (model 2). (4) 

iv) 

where i^^ '^'i is the sum over the bond (ij) connecting the vertices i 

and j, and i'^ "^2 is the sum over bonds i and j, which contain not 

only bonds in the chains but also virtual bonds that connect the center and 
the corners of the rigid junctions; S2 is defined on the vertices including 
the corners of the rigid junctions. 

The size of the rigid junctions can be characterized by the edge length 
R, which is fixed to 

R~ 0.1 (edge length of the rigid junctions) (5) 

in the simulations. The edge length i? = O.I is smaller than the mean 
bond length (~ 0.707); the bond length 0.707 corresponds to that in the 
equilibrium configuration of surfaces without the rigid junctions, where the 
relation 5'i/iV = 1.5 is satisfied. The rigid junctions in FiglUb) were drawn 
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to have a size R that is comparable to the mean bond length; the size R in 
FigHJb) was drawn many times larger than R — 0.1. As we will see later, 
the relation S'i/A^=1.5 is almost satisfied in model 2. 

The partition function Z of model 1 and model 2 is defined by 

Z = I l[dX,cxp[-S{X)], (6) 

i—l 

where S{X) is the Hamiltonian, which is given in Eq.© or in Eq.Q. j' 
denotes that the center of the meshwork is fixed in the integration. In model 
1, the dynamical variables are integrated over SiV-dimensional multiple 
integrations J' YliLi'i^i: while in model 2 they are integrated over 3Ns- 
dimensional integrations j' Jlilfi '^-^i ^^r the vertices and 6A^j-dimensional 
integrations /' Yii^Ji for the rigid junctions; 

J n^^'= (y J^^^^v [J Jl'^^T 

where N^Ns+Nj. 

We must emphasize that the definitions of the models in this paper are 
quite different from those of the conventional surface models including the 
compartmentalized surface models such as the one in |29) : The conventional 
surface models are defined on the triangulated surfaces, which always in- 
clude triangles (or plaquettes) even if the surface shape is maintained only 
by the skeletons. On the contrary, the meshwork in this paper has no pla- 
quettes except at the junctions in model 2 and is composed of the linear 
chains and the junctions only. 



3 Monte Carlo technique 

In model 1, the integration of the dynamical variable X is simulated by the 
random 3-D shift from X to X' ^X + SX, where 6X is randomly chosen 
in a small sphere. The new position X' is accepted with the probability 
Min[l, exp((55)], where (5S' = 5(new) — S'(old). The radius of the small sphere 
5X is fixed to certain constant value at the beginning of the simulations so 
that the acceptance rate is equal to about 50%. 

The vertices on the junctions share an energy which is different from 
that shared by the vertices on the chains in model 1. For this reason, 
we adopt an additional random shift for the vertices on the junctions. 
The first step is a simultaneous 3-D random translation of vertices on a 
junction, and the second step is a simultaneous 3-D random rotation of 
those vertices. Both of the shifts are done under about 50% acceptance 
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rate. Note that these additional shifts of X are not always necessary to 
simulate the integrations in model 1. 

The integration of X in model 2 is performed by a random 31? shift 
of the vertices on the chains and random 31? shifts of the rigid junctions, 
which respectively corresponds to the integrations in Eq. ([7]) . The 31? shifts 
of the rigid junctions are done by a random 3Z? translation and a random 
2>D rotation of the rigid plates. All of the shifts of X are done respectively 
under about 50% acceptance rate. 

A random number sequence called Mersenne Twister |32| is used in the 
simulations. 



4 Results of simulation 



(a) Collapsed (b) Smooth (c) Collapsed (d) Smooth 

6 = 2.92 6 = 2.93 6 = 2.59 6 = 2.61 

(model 1) (model 1) (model 2) (model 2) 



(e) The section (f) The section (g) The section (h) The section 



Fig. 3 (Color online) Snapshots of model 1 of size {N,Ns,Nj,L) = 
(36002,20250,2252,4) obtained at (a) 6 = 2.92 (collapsed phase), (b) 6 = 
2.93 (smooth phase), and snapshots of model 2 of size {N,Ns,Nj,L) = 
(22502,20250,2252,4) obtained at (c) 6 = 2.59 (collapsed phase), (b) 6 = 2.61 
(smooth phase). The mean square size , which is defined by Eq.®, is (a) 
X2~549, (b) X2~1563, (c) c^290, and (d) ~ 832. (e), (f), (g), and (h) 
are the meshwork sections. 

First we show that both meshwork models have a two-dimensional sur- 
face structure at least in the smooth phase. Figures [3]^ a) and O^b) are 
snapshots of model 1 surface of size {N, Ns,Nj,L) = (36002, 20250, 2252, 4) 
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obtained a.t b — 2.92 (collapsed phase) and b — 2.93 (smooth phase), and 
FigsEl^c) and [^d) are those of model 2 surface of size {N,Ns,Nj,L) — 
(22502,20250,2252,4) obtained at & = 2.59 (collapsed phase) and b = 
2.61 (smooth phase). The corresponding meshwork sections are shown in 
Figsl3l^e)-[3jh). The snapshots of model 1 were drawn in the same scale, 
which is slightly different from the scale for model 2. We see that the 
smooth state of model 1 in FiglSKf) contains the empty space inside the 
meshwork just like the conventional surface model '26]. However, the empty 
space is almost invisible in the crumpled state in FiglSjc). On the contrary, 
both of the smooth state and the collapsed state of model 2 in FigslS^g) 
and[3Kh) have empty spaces inside the meshworks. 



Si/N 
1.51 



1.49 



(a) 



model 1 



O :N=36002 
A:N=10242 
□ :N=2562 



Si/N 



1.51 



2.9 



1.5 

1.49 
(b) 



model 2 



O :N=22502 
A :N=6402 
□ :N=1602 



-o- 



2.6 



2.7 



Fig. 4 The Gaussian bond potential Si/N versus b of (a) model 1 and (b) model 
2. Si/N satisfies the predicted relation Si/A'^~1.5 in model 1, and it is slightly 
deviated from 1.5 in model 2. The solid lin es c onnecting the data were drawn by 
the multihistogram reweighting technique [34j . 



The relation Si/N = 3{N — 1)/2N ~ 1.5 is expected to be satisfied in 
model 1 because of the scale invariant property of the partition function 
|33j . On the other hand, the expected relation can slightly be violated in 
model 2. This is because of the finite size of the rigid junctions, although 
the scale invariant property is still valid in model 2. Figures lU^a) and|4l(b) 
show Si/N versus b of model 1 and model 2. We can see from the figures 
that the expectations are fulfilled. 

Therefore, we understand from the results shown in Figlll^a) that the 
MC simulations for model 1 were successfully performed. We consider that 
the MC simulations for model 2 are as well, because the simulation tech- 
nique for model 2 is almost identical to that for model 1. We find also from 
the results in FigsHJa) andlDJb) that the phase transition is not reflected 
in Si/N in contrast to the fluid surface model in Ref . [30j . where Si/N 
discontinuously changes at the transition point. 
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Fig. 5 The mean square size of versus h of (a) model 1 and (b) model 2. The 
curves were drawn by the multihistogram reweighting technique. 

The mean square size X"^ is defined by 

i i 

where X is the center of the meshwork. X"^ represents the distribution of 
vertices in R^. Therefore, the meshwork is expected to be characterized by 
large X"^ (small X^) at &^cxd(&-^0), where the meshwork is in a smooth 
(collapsed) state. 

Figures [5{a) and [H^b) show X^ versus h of model 1 and model 2. 
The solid curves drawn on the data were obtained by the multihistogram 
reweighting technique [34^. We find from FiglSl^a) that X^ changes almost 
discontinuously at an intermediate 6. Although the multihistogram curves 
at the transition point appear to be smooth, the discontinuous change of 
X"^ is apparent in the N = 19362 and N = 36002 surfaces. This indicates 
the existence of a discontinuous phase transition of surface-collapsing phe- 
nomenon between the smooth spherical phase and the collapsed phase. 
On the contrary, X^ of model 2 in FiglSjb) appears to vary continuously 
against b, and therefore, we can not confirm the discontinuous nature of 
the transition in model 2 only from X^ in Figl5l^b). 

In order to see the order of the collapsing transition in both models, we 
compute the variance Cx^ of X^ by 

Cx^.^j^{iX^~{X')f), (9) 

which can reflect how large the skeleton size fluctuates. If the models un- 
dergo a collapsing transition, we can see an anomalous peak in Cx^ at the 
transition point. We plot Cx^ versus b in FigslHI^a) and[6jb), which were 
obtained in model 1 and model 2, respectively. 
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The solid curves in FigsEl^a) and Et^b) were drawn by using the mul- 
tihistograni reweighting technique as weU as the corresponding curves 
in Figsin^a) and[5fb). We find in the curve of Cx'^ the expected anoma- 
lous peak, which apparently increases with increasing N . This indicates the 
existence of the collapsing transition in both models. 

The order of the collapsing transition can be confirmed from the scaling 
of the peak values C™^'^ such that 

C'^r oc {Nf' , (10) 

where (Ti is a critical exponent. We show in FiglHIc) the log-log plots of 
C^f^ against N , which were obtained from the curves in FigslHI^a) andlljb). 
The straight hues in Figlgl^c) were drawn by fitting the data to Ea.(fTO|. 
Thus, we have 

CTi = 1.75 ±0.05 (model 1), cti = 1.55 ± 0.09 (model 2). (11) 

The finite-size scaling (FSS) theory indicates that the peak values C^f^ 
should scale according to N'^{a=\) if the transition is of first order [35 j l36 1 
I37j . The exponents ui in both models are larger than 1, however, the ob- 
served scaling behavior in FigElJc) clearly reflects the anomalous behavior 
in Cx2. Therefore, the FSS analysis confirms that the model 2 undergoes a 
discontinuous collapsing transition, which was not always confirmed from 
the variation of X"^ against b in Figl^Kb). Moreover, the FSS analysis of 
model 1 is considered to be consistent to the first-order collapsing transi- 
tion indicated by the discontinuity of in FiglSl^a) . 

The discontinuous collapsing transition in both models can also be seen 
in the variations oi X^ . Figures [7{a)-[7I^i) show the variation oi X^ against 
MCS obtained at h which are close to the transition point of model 1. 
The surfaces are of size N = 10242, N = 19362, and N = 36002. We find 
that X'^ in Fig[7fa) remains in a lower value compared to that in FiglTjc). 
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Fig. 7 The variation of versus MCS of model 1 surface of size N = 10242 
at (a) & = 2.9, (b) b = 2.94, and (c) b = 2.98, and those of size iV = 19362 at (d) 
& = 2.92, (e) 6 = 2.93, and (f) 6 = 2.94, and those of size iV = 36002 at (g) 6 = 2.92, 
(h) 6 = 2.93, and (i) 6 = 2.94. Horizontal dashed hnes denote X^^^ and X'^^^ for 
computing the mean value of X^ in the collapsed phase and in the smooth phase, 
and and X^ax are shown in Table [1] The symbols smo and col in (b) denote 

the smooth phase and the collapsed phase, respectively. 



The large fluctuation of in Fig[7I^b) is consistent to the discontinuous 
transition between the collapsed phase and the smooth phase, which are 
respectively characterized by in FigslTlJa) and FiglTlJc). Contrary to 
of the N = 10242 surface, we find no such fluctuation in X^ on the 
iV = 19362 and = 36002 surfaces. After the surface is once trapped in the 
collapsed phase, it hardly changes to the smooth phase in those large sized 
surfaces. For this reason, a single variation of X^ in the collapsed phase and 
two variations of X'^ in the smooth phase are shown in Figsl7jd)-[7l^f) and 
[71Jg)-[7{i) . Horizontal dashed hnes denote X^-^^ and X^^^ for computing 
the mean value of X^ in the collapsed phase and in the smooth phase. The 
Hausdorff dimension can be extracted from these mean values of X^ and 
will be discussed below. 

The variations of X'^ of model 2 are shown in Figsl5]Ja)-[5^i) obtained 
at b which are close to the transition point of surfaces of size N = 6402, 
= 12102, and A^ = 22502. We find from the figures that the surfaces have 
two distinct states, which are respectively characterized by a large X^ and 
a small X'^ just the same as those in Figsl71[a)-[7Ui) of model 1. However, 
jumps of X^ can be seen in the variation in FigsO^e) and [Hh) even on 
such large surfaces. This indicates that the collapsing transition of model 
2 is relatively weak compared to that of model 1. Horizontal dashed lines 
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Fig. 8 The variation of versus MCS of model 2 surface of size A'' = 6402 at 
(a) b = 2.58, (b) 6 = 2.6, and (c) 6 = 2.62, and those of size iV = 12102 at (d) 
6 = 2.58, (e) 6 = 2.6, and (f) 6 = 2.62, and those of size iV = 22502 at (g) 6 = 2.59, 
(h) 6 = 2.6, and (i) 6 = 2.61. Horizontal dashed lines denote X^^^ and X^^ax for 
computing the mean value of X^ in the collapsed phase and in the smooth phase. 



denote X^^^ and X^ax f^'' computing the mean value of in the collapsed 
phase and in the smooth phase, just as in Figsl71[a)-[7Ji). 

The Hausdorff dimension H, which was already introduced in the final 
part of the Introduction, is obtained from the scaling relation 

X^ oc N^/" (12) 

by using the mean values obtained in the collapsed phase and in the 
smooth phase. 

As indicated by the horizontal dashed lines in Figsl7fa)-[7Ui) and FigslSfa)- 
[HUi), a lower bound X^f°^ and an upper bound X^™' ™ collapsed phase 
and those X^f™° and in the smooth phase can be assumed for com- 

puting the mean values of X'^. Table [T] shows the assumed values for X^^^ 
and ^max- Then, the mean values can be obtained from X'^ in the ranges 
^mfn' < < ^max and < X^ < The symbols 6(col) and 

6(smo) in Table [1] denote that the bending rigidity where the sequence of 

was obtained for computing the mean value of X^. The integers {£, m) 
are those introduced in Section[2]and characterize the size of the meshwork. 

FigureslHa) and[5](b) show the log- log plots of the mean values X^(smo) 
and X^(col) against N of model 1 and model 2, respectively. The error bars 
are the standard deviations. 
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Table 1 The lower bound and the upper bound Xn^ax' for computing the 

mean value X'^(col) in the collapsed phase, and those X^f™° and X^^° in the 
smooth phase. 6(col) and 6(smo) denote the bending rigidities where were 
obtained. 



model 


N 


{i,m) 


6(col) 


y-2 col 
^min 


y-2 col 


b(smo) 


V-2 smo 
^min 


Y'^ smo 
^max 


1 


5762 


(36,6) 


2.9 


80 


210 


2.98 


230 


360 


1 


10242 


(48,8) 


2.94 


110 


300 


2.98 


400 


600 


1 


19362 


(66,11) 


2.92 


200 


550 


2.93 


650 


1050 


1 


36002 


(90,15) 


2.92 


400 


730 


2.93 


1250 


1800 


2 


3602 


(36,6) 


2.6 


52 


123 


2.64 


140 


205 


2 


6402 


(48,8) 


2.58 


60 


175 


2.62 


210 


330 


2 


12102 


(66,11) 


2.58 


100 


300 


2.62 


400 


600 


2 


22502 


(90,15) 


2.59 


190 


500 


2.61 


680 


1000 




5000 10000 5000 10000 

(a) N (b) N 

Fig. 9 Log- log plots of X^(col) and X^(smo) against A'^ of (a) model 1 and (b) 
model 2. The straght lines were drawn by fitting the data to Ea. (|12p . The fitting 
in the collapsed phase in (b) was done by using the three largest data. 



The straight lines on the figures were obtained by fitting the data to 
Eq. (|12p . The fittings were done by using four data in the figures except the 
case of the collapsed phase of model 2 in Fig[9l[b). Thus, we have 

H''°^ = 2.68 ± 0.45, 77™^° = 2.15 ± 0.15 (model 1), 

77==°' = 2.31 ± 0.75, 2.20 ±0.15 (model 2). (13) 

We find from the results in Eq. (fT3l) that H^'^° in the smooth phase of 
both models are almost identical to the topological dimension H = 2, and 
that H"^"^ of model 1 in the collapsed phase is different from H^"'^° in the 
smooth phase, while H'^°^ and H^'^° of model 2 are almost identical. H"^"^ — 
2.31(75) of model 2 implies that the meshwork is relatively smooth even in 
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the collapsed phase, although the meshwork size discontinuously changes 
at the transition point. This is very similar to the case of the surface model 
with many holes in Ref.l38j- To the contrary, 7?^=°' = 2.68(45) of model 1 
implies that the meshwork is considerably collapsed in the collapsed phase. 
However, both H'^°^ remain in the physical bound, i.e., H'^°^<3. We should 
note that the results H'^°^ of both models are in sharp contrast to that of 
the compartmentalized model in Ref.[29], because the compartmentalized 
surface is completely collapsed in the collapsed phase. 



S2/N S2/N' 




Fig. 10 (a) The bending energy S2/N versus b of model 1 and (b) the bending 
energy S2/N' versus b of model 2. S2 is defined by Eq.© and Eq.®. iV and iV' 
are the total number of vertices where 5*2 is defined on the surfaces of model 1 
and model 2, respectively. 



We are interested also in the surface fluctuation phenomena; it is in- 
teresting to see whether or not the collapsing transition is accompanied by 
the surface fluctuation transition. The surface fluctuation phenomena can 
be characterized by the bending energy S2 in Eq.(l2l) and Eq.(j4|), because 
S2 is expected to be large (small) in the fluctuated (smooth) state in the 
meshwork model, just as the standard two-dimensional bending energy in 
the surface model. 

To see the surface fluctuations, we plot the bending energies iS'2/iV and 
S2/N' in Figs flOy a) and fTOT b). respectively. N and N' are the total number 
of vertices where the bending energy is defined in model 1 and model 2, 
respectively. N' in model 2 is given by N' = Ns + QNj — 12 = 30m{£ — m). 
Thus, we find that the variation of S2/N is discontinuous against b in 
model 1; the jump can be seen in S2/N on the surfaces of N = 19362 and 
N = 36002, although the gap is relatively small compared to the value of 
S2/N itself. From the discontinuous change in S2/N, we consider that the 



16 



surface fluctuation transition is of first order in model 1. On the contrary, 
S2/N' of model 2 in Fig fTOT b) appears to vary continuously against b. 




Fig. 11 The specific heat Csa for S2 versus b of (a) model 1 and (b) model 2, and 
(c) log-log plots of the peak values Cs^"" versus N{N') of model 1 (model 2). The 
error bars on the symbol denote the statistical errors, which were obtained by the 
binning analysis. Solid curves in (a) and (b) were drawn by the multihistogram 
reweighting technique. The straight lines in (c) were drawn by fitting the largest 
thr ee d ata (model 2) and the largest four data (model 1) to the expressions of 
Eq.(IT5}. 



The specific heat for 5*2 of model 1 is defined by 

Cs.=^{iS2~{S2)f), (14) 

which is the variance of ^2. of model 2 can be obtained by replacing N 
by N' in Eg. p^ . Just as Cx^ that reflects the collapsing transition, 
can also reflect the surface fluctuation transition through its anomalous be- 
havior. Figures fTlT a) and fTlT b) show versus b of model 1 and model 2, 
respectively. The curves on the figures were obtained also by the multihis- 
togram reweighting technique. We find an anomalous peak in of both 
models in Figs fTlT a) and [TlT b). The peak values increase with increasing 
N {N') and, therefore this indicates the existence of phase transition in 
both models. 

In order to see the order of the surface fluctuation transition, we show 
in FigdHc) the log-log plots of the peak values Cg'^'''^ against N{N'); N 
for model 1 and N' for model 2. Cg]^^ were obtained from Figs fTlT a) and 
ITlT b). The straight lines were drawn by fitting the data to 

eg'''' oc (Ny (model 1), Cg'''' oc {N' f (model 2), (15) 

where a2 is a critical exponent. The fittings were done by using the largest 
four data in model 1 and the largest three data in model 2 in FigfTlTcl. 
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Thus, we have 

(72 = 0.85 ±0.14 (modell), da = 0.51 ± 0.03 (model 2). (16) 

The exponent a2 of model 1 can be seen as (72 ~0.99 within the error and, 
hence it is almost equal to 1. Therefore, the scaling property of C^^^ seems 
consistent to the discontinuous transition indicated by the discontinuity of 
S2/N in Fig fTOT a). On the contrary, the exponent (72 of model 2 is obviously 
smaller than 1 and, therefore, this indicates that the surface fluctuation 
transition is considered to be of second order in model 2. 



Bx2 

0.6 
0.4 
0.2 




(a) 



(c) 



N=36002 



model 1 



model 1 



2.9 2.95 b 3 

0.6 
0.4 
0.2 




N=36002 



2.9 2.95 , 3 

b 



Bx2 

0.6 
0.4 
0.2 




(b) 



N=22502 



model 2 



2.55 2.6 ^2.65 2.7 




model 2 



N=22502 



2.55 2.6 ^2.65 2.7 



Fig. 12 The Binder quantity versus b of (a) model 1 and (b) model 2, and 
versus b of (c) model 1 and (d) model 2. 



The order of the transitions can also be characterized by the Binder 
quantities Bx^ and [39^, which are defined by 

<(^"-'^'»;>, B.. ^ 1 - - i^'":' . (17, 

If the transition is of first order, then we should have B = 2/3. In Figs fl2r a)- 
[T^ dl we plot Bx2 and Bg^ of both models. The solid lines in the figures 
were drawn by the multihistogram reweighting technique. 

We find in FigsfTWa) and [TWb) that both Bx^ and Bs2 have a peak 
j^max g^-j- ^-|-^^^ transition point and that B'^2^ is close to 2/3 at the transi- 
tion point of model 1, and also in Fig fT^ b) that B'^2^ is close to 2/3 at 
the transition point of model 2. These results on the order of the collapsing 



18 



transition in both models are consistent with the predictions by the discon- 
tinuity of in FiglUa) and by the finite size scahng analysis of C^f^ in 
Figslll^c). It is also seen that B^^^ of model 1 in FiglT2fc) is relatively close 
to 2/3 and that Bs^ of model 2 in FigflWd) is relatively smaller than 2/3, 
and therefore these results are also consistent with the results predicted by 
the observations in the FSS analyses of Cgj. 



5 Summary and Conclusion 

We have studied numerically two types of surface models defined on mesh- 
works, which are constructed as a sublattice in a triangulated surface and 
are composed of linear chains and junctions. It was found that both mod- 
els undergo a discontinuous collapsing transition between the smooth phase 
and the collapsed phase. The collapsed phases in both models are physi- 
cal in the sense that iJ < 3, where H is the Hausdorff dimension. With 
respect to surface fluctuations, the first model undergoes a discontinuous 
transition, while the second model a continuous transition. 

More precisely, two types of meshwork models of spherical topology were 
investigated by MC simulations for clarifying the phase structure, which is 
the order of the collapsing transition and the order of the surface fluctuation 
transition. Both models are deflned on meshworks, which are composed of 
linear chains and junctions; no two-dimensional surface is included in the 
meshwork except at the junctions. The flrst model, denoted by model 1, is 
characterized by elastic junctions, which are composed of vertices, bonds 
and triangles, and the shape of the elastic junctions is of hexagonal and 
of pentagonal. The Hamiltonian of model 1 contains the Gaussian bond 
potential, the one-dimensional bending energy, and the two-dimensional 
bending energy at the elastic junctions. The second model, denoted by 
model 2, is characterized by rigid junctions, which are hexagonal rigid plates 
and pentagonal ones. The Hamiltonian of model 2 contains the Gaussian 
bond potential and the one-dimensional bending energy. 

The bending rigidity bj at the elastic junctions was fixed to bj = 5 in 
model 1, and the edge length R of the rigid junctions was assumed to be 
R — Q.l in model 2. The compartment size was assumed to be i = 4, which 
is the total number of bonds in a chain between the junctions. Thus, the 
compartment size can be negligible compared to the surface size if N is 
sufficiently large. We used the lattices of size up to = 36002 in model 1 
and those up to TV = 22502 in model 2. 

We found that model 1 undergoes a first-order collapsing transition and 
a first-order surface fluctuation transition, and that model 2 undergoes 
a first-order collapsing transition and a second-order surface fluctuation 
transition. The smooth phase in both models is characterized by Hausdorff 
dimension ~ 2, while the collapsed phases are slightly different from 
each other between the two models. The Hausdorff dimension H"^"^ in the 
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collapsed phase at the transition point of model 1 is relatively larger than 
the topological dimension of the surface, but it remains in the physical 
bound, i.e., H <3. To the contrary, H'^°^ in the collapsed phase of model 
2 is almost identical to in the smooth phase. This implies that the 

collapsed meshwork of model 2 has a two-dimensional surface structure 
just as in the smooth phase although the meshwork size discontinuously 
changes at the transition point. 

Our results in this paper also indicate that the phase structure of the 
meshwork model is dependent on the elasticity at the junctions. In fact, 
the surface fluctuation transition is of first-order in model 1 while that is 
of second-order in model 2. 

Finally we comment on the relation between the models in this paper 
and the conventional fixed connectivity surface model [26| . We consider that 
the meshwork models in this paper are almost identical to the conventional 
surface models because of the following three reasons: Firstly, both the 
conventional surface model and the meshwork models have a discontinuous 
collapsing transition. Secondly, both the conventional model and model 
1 in this paper undergo a discontinuous transition of surface fluctuation. 
Thirdly, the collapsed phase of the meshwork models in this paper are 
physical, i.e., H <3, just as in the conventional model in [26], although no 
self-avoiding property is assumed in those models. 
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